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Abstract 

We study the irreversible dynamics of nonlinear, nonintegrable Hamiltonian oscil- 
lator chains approaching their statistical asymptotic states. In systems constrained 
by more than one conserved quantity, the partitioning of the conserved quantities 
leads naturally to localized and coherent structures. If the phase space is compact, 
the final equilibrium state is governed by entropy maximization and the coherent 
structures are stable lumps. In systems where the phase space is not compact, the 
coherent structures can be collapses represented in phase space by a heteroclinic 
connection of some unstable saddle to infinity. 
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1 Introduction 



The formation of coherent structures by a continuing self focusing process is a 
widespread phenomenon in dispersive nonlinear wave systems. As a result of 
this process, high peaks of some physical field emerge from a low amplitude 
noisy background. For optical waves in Kerr-nonlinear media this results in a 
self-enhanced increase of light intensity in a small sector of a laser beam while 
the intensity in the neighborhood of this bright spot decreases [1]. Related 
phenomena occur in diverse systems from hydrodynamics to collapsing Lang- 
muir waves in a plasma [2] and in numerical algorithms for partial differential 
equations [3], [4] . 
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A common feature of the wave dynamics of these systems is the comparable 
strength of the dispersion and the nonlinearity. But, self focusing phenomena 
are radically different in integrable and nonintegrable systems. In integrable 
systems [5], the peaks appear and disappear in a quasiperiodic manner re- 
flecting the phase space structure of nested tori. This behavior is usually 
encountered on time scales that are short enough so that generic noninte- 
grable contributions to the dynamics may be neglected. Significant changes 
occur on time scales where the nonintegrability is relevant [6], [7]. The solu- 
tion's shape becomes more irregular [8], [9], [10], [11] and the periodic breathing 
of the peaks turns into a more persistent state. These peaks can merge into 
stronger ones while radiating low-amplitude waves. The irreversible charac- 
ter of the system becomes apparent and its behavior is driven by statistical 
mechanics [12], [13], [14] as the solution's trajectory tries to explore more of 
the available phase space. As it does this, it must, at the same time, respect 
the preservation of conserved quantities. The phase space shell of constant 
conserved quantities determines the system's most favorable macrostate and 
the resulting dynamics can lead to a local gathering of the amplitude while 
the system explores the phase space shell. We will see in situations with more 
than one conserved quantity that the spatial structure of the state which is 
statistically preferred now contains coherent peaks. If the phase space is com- 
pact, these peaks can be local , time independent and stable. If the phase 
space is not compact, the peaks can be collapsing filaments which produce 
singularities. 

In this paper we study three systems which typify the focusing behavior ob- 
served in nonintegrable Hamiltonian systems with more than one conserved 
quantity, namely the Landau-Lifshitz equation for a classical Heisenberg spin 
chain, various versions of the discrete nonlinear Schrodinger equation, and a 
leapfrog discretization of the Korteweg-de Vries equation. Various modifica- 
tions of these systems help to identify the role of the conserved quantities 
in the formation of coherent structures. The spatial discreteness of all these 
systems avoids fluctuations on infinitesimal scales. 

The first case study investigates the Landau-Lifshitz equation for the classical 
spin chain in one spatial dimension. We study the long time behavior of those 
configurations in which most spins are close to the north pole. The two con- 
served quantities are the energy and the magnetic moment. Almost constant 
states undergo a series of modulational instabilities and the system begins 
to oscillate as if it were integrable. Nonintegrability, however, leads to non 
recurrence as localized peaks appear which merge from time to time leading 
to even larger peaks and radiating some energy. The phase space is compact, 
and one can compute to good accuracy the thermodynamic potentials of the 
system by separating the low-amplitude spin waves and the strongly nonlinear 
components. Depending on the initial value of the energy and the magnetiza- 
tion, entropy maximization leads to a state where some of the magnetization 
must be put into local structures bounded by domain walls for which the spin 
of each contained lattice point is close to the south pole. Numerical simula- 
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tions very clearly support our simple analytical predictions which are based 
on thermodynamic considerations. 

A close relative of the Heisenberg spin chain is found by taking the small 
amplitude limit where all spins are close to the north pole. Deviations are de- 
scribed by the focusing nonlinear Schrodinger equation. A similar dynamics is 
observed. In this study, we can also investigate the effects of exact integrability 
by using the Ablowitz-Ladik algorithm [15]. In these simulations, we find no 
irreversible behavior, no peak fusion, no relaxation, only quasi-periodic behav- 
ior. Likewise, we get qualitatively different results if we take a model which 
breaks the rotational symmetry and because of this the particle number is 
no longer conserved. As a consequence it is not necessary for the system to 
develop coherent structures in order to maximize its entropy as it no longer 
has to be concerned about the second constant of motion when its trajectory 
explores the accessible phase space. 

The last case study concerns the leapfrog algorithm for the numerical inte- 
gration of the Korteweg-deVries equation. It was observed in previous works 
[3], [4] that the leapfrog algorithm for the Korteweg-deVries equation always 
develops singularities. After (usually) a very long time, the amplitudes in some 
local neighborhood rapidly diverges. We demonstrate that this collapse again 
takes on an organized coherent form. 

How and why does the system develop such local objects? The reason is again 
statistical. At an early stage, one can again observe the gathering of one of 
the conserved quantities in coherent structures. The system's phase space is 
not compact, however, so that a strong nonequilibrium process prevails finally. 
We find a rapidly growing localized 'monster'-solution (so called because of 
its likeness in shape to the Loch Ness monster) that has a canonical structure. 
This solution can originate from coherent structures or, most frequently, from 
a long wave instability of the low-amplitude noisy background. We discuss 
the similarity of this process to the collapse behavior of the two dimensional 
focusing nonlinear Schrodinger equation, where the conservation laws necessi- 
tate a net inverse particle flux to small wavenumbers. 

This paper is arranged as follows: In section 2, we present the Landau-Lifshitz 
equation, the discrete nonlinear Schrodinger equation and the leapfrog dis- 
cretization of the Korteweg - de Vries equation and discuss some of their 
properties. In section 3 we present numerical studies of these equations. In par- 
ticular, we study thermodynamic quantities during the focusing process. Their 
significance is also demonstrated by discussing modified equations with either 
more or with fewer conserved quantities as well as an equation of defocusing 
type. In section 4 we will give a statistical interpretation of the numerical find- 
ings. By computing the thermodynamic potentials of the spin-system, we find 
the connections between global features of the pattern and the conserved quan- 
tities. Macroscopic properties of the final state are computed. The discretized 
KdV equation has no such state of thermal equilibrium. We identify the rapid 
divergence of the amplitudes with the exploration of the non-compact phase 
space shell and we suggest that there is much similarity between this behav- 
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ior and the condensation and collapse behavior seen in the focusing nonlinear 
Schrodinger equation. 

Figures of similar contents are grouped together and their order sometimes 
deviates from the sequence of their references in the text. 



2 Nonintegrable systems with constraints 

2.1 Time- continuous systems 



2.1.1 The Landau- Lifshitz equation 

The anisotropic Heisenberg spin chain is particularly suitable for the study 
of self-focusing phenomena. This system contains the generic properties of 
equations of nonlinear Schrodinger type that lead to self-focusing and it is 
easy to investigate from the statistical point of view. The Landau-Lifshitz 
equation [16] 

S n = S„ x (J(S„_i + S n+ i) + S nz e z ) (1) 



is a classical approximation of the dynamics of magnetic moments S n 

n. S n is perpendicular to S n . Therefore the 
moduli of the spin vectors are conserved and one may set |S n | = 1. The 
phase space of a chain of N spins is a product of iV such spheres. S z the 
component of S along the rotational symmetry axis. The northern and south- 
ern hemispheres are equivalent since (1) is invariant under the transformation 

(S X ni Sy n , S Z n) ^ ( S xn , Sy n , S zn ). 

There are two trivial homogeneous equilibrium states where all the spins point 
either to the north pole S z — 1 or to the south pole S z = —1. Throughout 
this paper we only consider solutions where most of the spins are close to the 
north pole. 



2.1.2 The discrete nonlinear Schrodinger equation 

The long-wavelength dynamics of spins which deviate slightly from the north 
pole S z — 1 is given by the focusing discrete nonlinear Schrodinger (DNLS) 
equation 

i(j) n = J(0 n+ i + 4> n -i - 2<p n ) + 2\<f) n \ 2 <f) n (2) 



for small values of the complex amplitude = (S x + iS y )/(l + S z ). The 
spin chain may thus be regarded as a DNLS which is modified by higher order 
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terms. The north pole is corresponds to = while the south pole corresponds 
to an infinite amplitude. 



2.1.3 Integrals of motion 

The spin chain and the DNLS equation each have two conserved quantities: 

(1) The Hamiltonian of the DNLS equation Ti = J2 n J (14 > n4>n ~ 0n0n+i ~~ 
0*0 n+ i) — \4> n \ A is again obtained as the lowest order of the Hamiltonian 
of the spin chain H = Hj + H a = En J(l ~ S n S n+ i) + (1 - S 2 zn )/2. The 
first contribution is a Heisenberg exchange coupling which is minimal for 
homogeneous solutions. The second part is an anisotropic energy that has 
minima at the poles S z — ±1 and is maximal at the equator S z = 0. The 
stationary spin-up or spin-down solutions are the absolute energy min- 
ima. Fluctuations about these ground states give energy contributions 
per lattice site of the order of \S X + iS y \ 2 . Similarly, small fluctuations 
near by the equilibrium state <fi = of the DNLS contribute a coupling 
energy proportional to |0| 2 . In contrast to the south pole state of the 
spin system, the energy of a solution <fi — > oo in the DNLS goes to minus 
infinity as — 10| 4 . 

(2) The second conserved quantity of each system, the total magnetization 
■M — Szn of the spin chain and the modulus-square norm ('particle 
number') Y. n \4>n\ 2 of the DNLS, are related to the system's rotational 
symmetry. The superposition of the Hamiltonian and this integral of 
motion yields a Hamiltonian in a rotating frame system. The negative 
magnetization N — M. = — S zn ) ~ \S X n + iS yn \ 2 /2 corresponds to 
the particle number of the DNLS in the lowest order in amplitude. This 
'particle number' of the spin chain is zero for the north pole solution and 
it is two per lattice site for the south pole solution. In the DNLS, the 
particle number diverges for the state |0| — > oo. 

In order to contrast the generic behavior of such systems with those of (a) inte- 
grable systems and (b) systems not constrained by a second conservation law, 
we also consider two modified equations of motion. The first is the integrable 
Ablowitz-Ladik discretization of the one-dimensional nonlinear Schrodinger 
equation [15]). The second is an equation where the second integral is de- 
stroyed by a symmetry breaking field. 

Low energetic solutions just above the ground states can be characterized by 
the ratio of the two integrals of motion, i.e. the energy per particle. This reveals 
a major difference between the spin chain and the DNLS. For fluctuations near 
the north pole or near = 0, the particle number is of the order of the energy 
so that this ratio is of order one for both systems. Spin-fluctuations near the 
south pole have the same energy but the second conserved quantity ('particle 
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number') is much higher so that the energy per particle is proportional to 
\S X + iS y \ 2 and much less than 1. Thus the spin chain has two states with 
low energies per lattice site, one (S z ~ 1) with a higher energy per particle 
and one (S z ~ —1) with a low positive energy per particle. This well-defined 
condensate state of low energy and high particle density is a mayor advantage 
of the spin chain. 

In contrast, for infinitely high amplitude solutions of the DNLS that corre- 
spond to S z ~ — 1 solutions of the spin chain both the energy and the particle 
density go to infinity. The energy per particle diverges proportional to — 10| 2 . 



2.2 Time-discretized equation of motion 

2.2.1 The leapfrog-discretization of the Korteweg-de Vries equation 
The system of finite difference equations 

u m +i(n) = v m (n) 

+r(u m (n + 2) - 2u m {n + 1) + 2u m {n - 1) - u m {n - 2) 
-2(u m (n + 1) + u m (n) + u m (n - l))(u m (n + 1) - u m (n - 1))) 

v m+ i(n) = u m (n) 

is a leapfrog-type discretization in space and time of the completely integrable 
Korteweg-de Vries (KdV) equation u = u xxx — 6uu x for the real amplitude 
u(x,t). The leapfrog discretization is characterized by a central difference 
du/dt — > (u m+ i — u m _i). The factor r of the spatial derivative is the time 
step size. The precursor w m _i is identified with the additional variable v rn on 
the right side of the first equation. This scheme allows the simulation of the 
partial differential equation avoiding the amplitude dissipation that occurs 
in methods with numerical viscosity. The term u m {n + 2) — 2u m (n + 1) + 
2u m {n — 1) — u m {n — 2) is the standard discretization of u xxx . The discretiza- 
tion of uu x was first suggested by Zabusky and Kruskal [17]. It involves a 
central difference in space u m (n+ 1) — u m {n — 1) and replaces u by the average 

+ 1) + u m (n) + u m (n — l))/3. 
This discretization suppresses fast-acting nonlinear instabilities. Discretiza- 
tions that do not retain some of the original conservation laws lead to fast 
acting instabilities, since single modes diverge rapidly. For instance, the mode 
with the wavenumber k = 27r/3 is driven by the nonlinear part of conventional 
discretizations of the KdV equation. In contrast, this mode is an exact solution 
of (3). Linear instabilities of the zero-solution can be avoided by a sufficiently 
small step-size r < 2/(3\/3). 
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2.2.2 Integrals of motion 

The special feature of the spatial discretization (3) is that it preserves some 
of the original conserved quantities: 

(1) < uv >= J2n u m{n)v m (n) = const, corresponds to the conserved quantity 
/ u 2 dx ('energy' for shallow water waves) in the original KdV equation. 
The modulus-square norm < u 2 + v 2 >= J2n u m( n ) 2 + v m (n) 2 is not 
conserved. 

(2) < U >= E„tt2m(n) = En^2m+l(n) and < V >= En^2m(n) = En u 2m+lW 

correspond to / udx ('mass' for shallow water waves). 



3 Numerical studies 

We examine the formation of coherent structures numerically in various ver- 
sions of the spin chain and the DNLS equation as well as the leapfrog inte- 
gration scheme for the KdV equation. A typical scenario for the spin chain 
suggests that the final state mainly depends on the amount of the two con- 
served quantities provided by the initial conditions. Simulations of various 
modifications of the DNLS with either more of less integrals of motion clarify 
some more general conditions for this behavior. The simulations of the differ- 
ential equations apply an Adams routine to a chain of 512 (and occasionally 
4096) oscillators with periodic boundary conditions. 

3. 1 Dynamics of the spin chain 

3.1.1 Benjamin-Feir instability and reversible dynamics 

Plane wave solutions of the nonlinear Schrodinger equation are Benjamin-Feir 
unstable so that self-focusing is initiated by long wavelength modulations. 
Similarly spin-wave solutions of the Heisenberg spin chain are unstable under 
long wavelength perturbations. For instance, the homogeneously magnetized 
solution S n = S (Fig. la) that precesses about the symmetry axis e z with the 
frequency uj — S z is most unstable under perturbations with the wavenumber 

k = vA - si/ j. 

As a result of this instability, small perturbations lead to spatially periodic 
humps of spins approaching the equator while most of the spins come closer to 
the north pole (Fig. lb). The trajectory is close to a homoclinic orbit so that the 
solution returns to the almost homogeneous state after reaching the maximum 
amplitude. Fig. 2a shows the profile of the z- magnetization in space and time 
for this solution which is almost periodic in space and time. Fig.3a shows the 
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Fig. 1. Sketch of the spins for (a) the spatially homogeneous solution, (b) the spa- 
tially periodic solution resulting from a Benjamin-Feir instability, (c) the humps 
moving towards each other merging into compound peaks of spins pointing down 
(d). 

fa) fb) 




Fig. 2. S zn as a function of the lattice site n and time, (a) portrays the sector (a) 
of Fig. 3a, (b) is the sector (/?); (c) is at (7) of Fig. 3b. 

pattern of peaks (the sites of the spins that differ most from the north pole) as 
a function of time; Fig. 2a is related to box (a) of Fig. 3a. Equidistant humps 
emerge and disappear periodically in time for t < 200. 
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Fig. 3. Integration of 512 spins (J = 0.4) with periodic boundary conditions; lat- 
tice sites where the spins deviate significantly from the north pole (S z < 0.8) 
are marked with dots. Integration over (a) 2000 and (b) 200,000 time steps with 
weakly perturbed homogeneous initial conditions (S z ps y/0.84) (the profile of the 
z-magnetization at (a, 0, 7) are given in Fig.2); noisy short wave (k = ir) inital con- 
ditions with a xenocryst of spins deviating strongly from the north pole over 2000 
time steps (c); defocusing equation (negative sign of the potential), weak coupling 
(J = 0.1) and noisy short wave (k = n) initial conditions over 20,000 time steps 
(d). 

3.1.2 Merging of peaks and irreversible dynamics 

The spatially periodic solution arising from the Benjamin-Feir instability is 
itself phase- unstable. As a result, the periodic pattern with the wavelength of 
the initial periodic mode becomes modulated on an even larger length scale 
so that the the gaps between the initially equidistant humps start to vary 
(Fig.lc). Fig.3a shows tiny variations of the distances between the humps as 
the humps start to move at t ~ 60 . 

The most important phenomenon following the phase instabilities is the for- 
mation of coherent structures through mergings of peaks. Neighboring humps 
approaching each other finally merge into single peaks radiating small fluc- 
tuations. Fig. 2b shows the profile of the magnetization during the fusion of 
humps of box (f3) in Fig. 3a. The original periodic solution is smooth, but the 
compound peak resulting from the merging has an irregular shape involving 
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huge gradients both in space and in time. Its amplitude oscillates irregularly 
in time, but unlike the original periodic solution, it does not vanish any more 
completely. Even those humps that are situated remotely from the first merg- 
ing processes become more persistent in time immediately so that they are 
traced by continuous lines in Fig. 3a. 

Subsequently, more humps fusing into compound peaks increase the average 
distance between neighboring peaks. The resulting compound peaks again 
merge with primary humps and with other compound peaks forming even 
stronger peaks (see point (7) in Fig. 3b with the magnetization profile of 
Fig.2c). 

3.1.3 Final equilibrium state 

The increasingly high amplitude of the compound-peaks enables some spins 
to overcome the energy barrier of the equator and to flip to the southern 
hemisphere (Fig. Id). (Fig.3b) shows that after 2 x 10 5 time steps all peaks 
have merged into six down-magnetized domains that each consist of two or 
three lattice points. These peaks with S z ~ — 1 are embedded in a disordered 
state where the spins deviate only slightly from the north pole S z = 1. The final 
state is a two domain pattern where most of the spins are accumulated in huge 
up-magnetized domains while a few spins condense to small down-magnetized 
domains. The domain of the spin-wave fluctuations near to the north pole in 
the final state remains persistent even if the spin-down xenocrysts are removed 
artificially by flipping the down-spins up as S z — > \S Z \. 

3.1.4 Transfer of energy 

The transition from the almost regular dynamics to the irreversible process 
during the merging of peaks is reflected in the share of the total energy of 
the two parts of the Hamiltonian. The coupling energy Tij = XX 1 — S n S n +i) 
results from spatial inhomogeneities within each of the two domains and from 
the domain walls between them. The anisotropic energy 7i a = |XXl — S^ z ) 
depends on the distance of the spins from the poles. Fig. 4a shows the transfer 
of energy between the two parts of the Hamiltonian: The initial state contains 
no coupling energy and the anisotropic term is the only contribution. Some 
of this energy flows to the coupling part during the formation of the spatially 
periodic pattern. This process is reversed while the system approaches the 
homogeneous state again so that energy is exchanged periodically between H a 
and 7ij (see the periodic behavior for short times in Fig.4a). 
The transfer of energy from the anisotropy to the coupling becomes irreversible 
when the humps fuse into compound peaks. The share of coupling energy 
increases even more when compound peaks merge and an increasing number of 
spins overcomes the equator and settles down near to the south pole. In phase 
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Fig. 4. (a): The evolution of the anisotropic energy E a and the coupling energy 
Ej over 2000 time steps for 4096 spins with the initial conditions of Fig. 3. (b): 
Numerical (points) and thermodynamic (line, see section 4.4) results of E a /Ej after 
long integration times. 

space, this process is related to Arnold diffusion. The trajectory disappears 
from the initial critical torus and explores more and more of the phase space. 
Finally, the system reaches an equilibrium state where the proportion of 7i a 
and Tij saturates (Fig.4b). The bulk contribution to the coupling energy is due 
to the inhomogeneity within the north-pole domains and not the contribution 
from the domain walls. 

3.2 Modified initial conditions 

The systems energy and magnetization per spin given by the initial conditions 
determines the number of spins that point down after a long time. We study 
the final state for initial conditions with varying energies and with a modified 
magnetization profile. This gives strong numerical indications that the two 
integrals of motion (energy and the magnetization) are the key quantities 
that influence the final state. 

3.2.1 Varying energies 

Spin-wave like initial conditions S xn + iS yn = ^Jl — S% exp(ikn) with a given 
amplitude provide energies E = N(l — 5f)(| + J(l — cosk)) depending on 
the wave number k while the magnetization M = NS Z is A;- independent. The 
initial condition k = in the simulation described in Fig. 2 (a), (b) corresponds 
to the minimal energy E — y(l — Si)/ 2 that is possible for a given magne- 
tization. The maximum energy E = N(l — Sf)(| + 2 J) corresponds to an 
excitation at the boundary of the Brillouin zone k = it and any energy be- 
tween these values can be obtained by a suitable spin wave. 
We find that some of the spins condense near to the south pole eventually only 
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if the energy is below a certain threshold E eq within this range. Fig. 5 shows 
the magnetization of the south-pole condensate as a function of the systems 
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Fig. 5. Numerical (dots) and analytical (line) results for the total magnetization 
J2n(s nz <0) \Snz\ of the spins in the southern hemisphere as a function of the to- 
tal energy E a + Ej per spin for the system of Fig. la. The initial conditions are 
S z = \/0.84, S xn + iSy n = 0.4exp(i/cn) with various wavenumbers k; waves with 
wavenumbers from k = to k = it provide different amounts of coupling energy Ej 
initially. The spin chain is integrated over 100000 time steps. For energies per spin 
below the threshold (equation (18)), some of the spins condense near to the south 
pole. 

energy: 

(i) Oversaturated phase: Below the energy threshold E eq , the number of 
down-spins is proportional to E eq — E. The behavior in this is very sim- 
ilar to the scenario described before. The spatially homogeneous initial 
condition in the above simulations just leads to the highest possible pro- 
portion of the south-pole condensate. 

(ii) Overheated phase: For high energies E > E eq , no spins are flipped 
down so that this magnetization is zero. There is no south-pole condensate 
beyond this threshold; all spins end up in small fluctuations near the 
north pole. 

3.2.2 Modified magnetization profile 

We have seen that only long- wavelength fluctuations create peaks. In contrast, 
high-energetic initial conditions with short wavelengths melt away such peaks 
of spins deviating significantly from the north pole. This occurs for an initial 
condition of a small-amplitude k = ir spin-wave (corresponding to the max- 
imum energy in Fig. 5) where the spins within a small domain are flipped to 
the southern hemisphere. Fig.3c shows the destruction of such a domain in a 
bath of k = n waves. The system ends up in an irregular state where all spins 
are near to the north pole. 
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3.3 Modified equations of motion 



The scenario we have described is widespread in dynamical systems and not a 
specific feature of the Landau-Lifshitz equation. A comparison of this scenario 
with self-focusing in related systems indicates that the main conditions for the 
emergence of coherent structures are 

(i) the low-amplitude dynamics is governed by an NLS-type of equation, 

(ii) the system is nonintegrable, 

(iii) there are two integrals of motion. 

The first of these points basically characterizes the dynamics of nonlinear dis- 
persive systems on long scales. The focusing DNLS is the system most closely 
related to the spin chain. Also the case of a defocusing nonlinearity will be 
considered. The importance of nonintegrability will be shown by the compar- 
ison to the integrable Ablowitz-Ladik discretization of the NLS equation. On 
the other side, we will study a system with broken rotational symmetry that 
only conserves the Hamiltonian. 

3.3.1 Defocusing equation 

The formation of coherent structures in discrete nonintegrable systems is not 
an exclusive property of 'focusing' types of discrete NLS or Landau-Lifshitz 
equations. The Landau-Lifshitz equation S n = S n x ((J(S n _i + S n+ i) — S nz e z ) 
with a negative ('easy-plane') anisotropy corresponds to the defocusing dis- 
crete NLS equation. For weak coupling constants (J = 0.1), short wavelength 
(k = 7r) initial conditions produce coherent structures with spins condens- 
ing in the equator region where the anisotropic energy has now its minimum. 
Fig. 3d shows this weak focusing process for the spin chain. 

3.3.2 Discrete NLS equation 

The focusing nonintegrable DNLS i<j) n = J(0 n+ i + 4> n -i — 20 n ) + 2|0„| 2 n 
has the properties (i)-(iii) just like the Heisenberg spin chain. Fig.6a shows 
the spatiotemporal pattern for the DNLS of lattice sites with high amplitudes 
that is very similar to the one described in section 3.1 (Fig. 2a). Again, an un- 
stable periodic pattern emerges from a phase-instability of the homogeneous 
state that is unstable itself. 

The first and the second phase instability are well-known as direct conse- 
quences of (i) and (ii). The second phase instability has been studied in detail 
in the context of NLS equations. In phase space, it is related to degenerate tori 
with less than the maximum dimension. Such critical tori exist in integrable 
as well as in nonintegrable systems; they may be stable or unstable. The sta- 
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bility of these tori is related to double points in the spectral transform [10]. 
The spectrum of the Lax-operators has been analyzed both for an integrable 
and a nonintegrable version of the discrete NLS equation. 
The fusions of humps lead to peaks with high amplitudes (Fig.6a) and fi- 
nally high-amplitude xenocrysts emerge from a low-amplitude turbulent back- 
ground. The whole process is very similar to the one of the Landau-Lifshitz 
equation (Fig. 3a). The correspondence persists even in a domain where the 



fa) fb) 




□ sun f man □ sun t man 



Fig. 6. Integration of the various versions of the discrete nonlinear Schrodinger 
equation (J = 0.4) with 512 lattice sites and periodic boundary conditions. The 
inital conditions are (ft n = 0.2 plus noise, lattice sites with \cf>\ > 0.25 are marked 
with dots, (a): Integration of the DNLS over 2000 time steps; (b): Integration of 
the of the integrable version of the DNLS over 2000 time steps, (c): DNLS with a 
contribution —2<p n and a symmetry breaking field 0.2<j)' n over 10,000 time steps; (d): 
DNLS with a symmetry breaking contribution 0.02^^ over 10,000 time steps. 

additional nonlinear terms in the Landau-Lifshitz equation are not small. The 
DNLS-peaks may have different heights while the spin-peaks are always south- 
pole states. 

The radiation of low-amplitude fluctuations during the merging of peaks is 
related to homoclinic chaos following the break up of Kolmogorov-Arnold- 
Moser tori. By discussing the Melnikov- function of the critical tori, [11] have 
detected homoclinic crossings in the nonintegrable NLS. The horseshoe of the 
homoclinic crossings creates the disorder following the fusion of two peaks. 
A phenomenon similar to the merging of peaks was found in a continuous non- 
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integrable NLS equation [7]. Unlike solitons in integrable systems, collisions of 
solitary solutions of nonintegrable equations lead to a transfer of power from 
the weaker soliton to the stronger one while low amplitude waves are radi- 
ated. The resulting two-component solution contains a decreasing number of 
growing solitons immersed in a sea of weakly turbulent waves. In continuous 
systems, energy is drained by infinitesimal scales while the spatial discretiza- 
tion defines a minimal length scale. 

3.3.3 Integrable discrete NLS equation 

Homoclinic chaos as a source of radiation is absent in integrable systems. 
The comparison with the integrable discrete NLS i<f) n = J(0 n+ i + </> n _i — 
20 n ) + |0n| 2 (0«-i + <f> n +i) shows this implication of the nonintegrability (ii). 
The integrable NLS equation exhibits the primary phase instability, but not 
the fusing of neighboring peaks. While the dynamics is similar to the nonin- 
tegrable system initially, the peaks do not merge (Fig. 6b). Consequently, no 
coherent structures evolve and the system does not settle down in a disordered 
equilibrium state. The quasiperiodic appearance of humps with relatively low 
amplitudes reflects the phase space structure of nested tori. 

3.3.4 Particle nonconserving equation of motion 

While the irreversible focusing process is a consequence of nonintegrability, 
it also depends on the existence of some remaining integrals of motion. The 
generation of coherent structures is very sensitive to perturbations that de- 
stroy one of the remaining integrals. The property (iii) may be changed by 
breaking the rotational symmetry with the contribution eRe(<p n ) in the NLS 
equation. The equation is still of Hamiltonian type, but the modulus-square 
norm X^|^n| 2 is not conserved. Additional contributions ~ ^J2\(pn\ 2 to the 
Hamiltonian and the corresponding term ~ uocf) in the equation of motion 
i(j) n = J((f) n+ i + (p n -i — 20„) + u(j) n + eRe((j) n ) + 2|0 n | 2 n are now relevant for 
the dynamics (in the symmetric case, this term just describes the same dy- 
namics in different rotating frame systems; in the symmetry broken system, 
the external field e is stationary in the system that rotates with the frequency 
a;). Depending on the sign of u, two different scenarios are observed: 

uj < 0: The onset of self-focusing for small times is similar to the symmet- 
ric case. However, the peaks emerging from the fusing process disintegrate 
eventually into small amplitude fluctuations (Fig. 6c). 

uj > 0: The onset of the focusing process is again similar to the one with 
particle conservation, but after about 5000 time steps growing amplitude 
fluctuations lead to a disordered state. Unlike the rotationally symmetric 
system, high particle densities are not confined to small islands in a sea 
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of low particle density fluctuations. The nonconservation of the particle 
number leads to high (but finite) particle density fluctuations everywhere 
(Fig.6d). 

In the spin chain, similar effects can be reached with an external magnetic 
field that is perpendicular to the anisotropy axis e z and an additional z-field. 
The Hamiltonian now contains the additional Zeeman terms eS x + uS z . Due 
to the broken rotational symmetry the total magnetization is no longer an 
integral of motion. 

3.4 The leapfrog- discretization of the Korteweg-deVries equation 

Iterations of the leapfrog-discretization of the Korteweg-deVries equation ex- 
hibit a scenario of merging peaks that is very similar to the one found in NLS 
or spin equations. However, the leapfrog system undergoes a rapid unbounded 
growth similar to the blow-up in twodimensional NLS systems. While the 
Heisenberg spin chain has a well-defined equilibrium, the leapfrog system al- 
lows us to study the conditions for blow-ups in constrained systems. We study 
this phenomenon for two initial conditions, the k = 2n/3 mode with a strong 
correlation of u and v, and for white noise with no correlation of u and v. The 
system consists of 1020 lattice sites with periodic boundary conditions. 

3.4-1 Correlated initial conditions 

The monochromatic wave with the wavenumber k = 2ir/3 as initial condition 
yields an exact but phase-unstable [18] solution of the leapfrog-iteration (3a) 
for a sufficiently small step size. This mode is particularly relevant for a sta- 
bility analysis since it is the fastest growing mode for a step size r > 2/(3\/3)- 
Setting Ui(n) = Vi(n) provides the maximum correlation of u and v. The 
pattern of peaks Fig. 7a (lattice sites with high u(n) 2 + v{n) 2 ) emerges in a 
manner similar to the spin- and NLS-systems (Fig. 3 and Fig. 6): 

(i) Regular behavior: The initial low-amplitude wave is below the thresh- 
old to be traced in Fig. 7 for m < 400. 

(ii) Merging humps: A phase instability of the initial k = 27r/3 wave leads 
to a modulational pattern with a wavelength of about 20 lattice sites 
that reaches the threshold atmw 400 so that a spatially periodic pattern 
emerges for 400 < m < 600. These periodic humps are wave packets of the 
initial short wave moving towards higher n. Similar to the spin- and NLS- 
systems, this pattern itself is slowly modulated. The humps approach 
each other and merge so that a decreasing number of peaks of increasing 
intensity survive. Solitary solutions that are high and fast sweep away 
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Fig. 7. Iteration of the leapfrog-discretized KdV-equation (3) for a k = 2ir/3 wave 
(a) and for white noise initial conditions (b) with periodic boundary conditions. 
The amplitude of u(n) = v(n) is 0.1 initially in (a). Locations with high amplitudes 
(the sum of u{n) 2 over five subsequent steps is greater than 0.1) are marked with a 
dot. In (b), the threshold of \/u 2 + v 2 is is 0.045 while the noise level ~ 0.01. 

slower ones. The peaks speed and amplitude increase while the width 
decreases during this process. They accumulate high amounts of the con- 
served quantity < uv > just like the spin-down domains gather magne- 
tization. Fig. 8 shows the cumulated conserved quantity Y%=i u m,(l) v m,(l) 




Fig. 8. Cumulated energy J2f=i Um(l)v m (l) for the simulation of Fig.7 and Fig.9 as 
a function of the lattice site n at the time steps m = 1 and m = 2500 

as a function of n (for n — N, it is conserved) at the beginning (m = 1) 
and at m — 2500 when the solitary waves have developed. While the con- 
served correlation J2u(n)v(n) is equally distributed in space initially, the 
formation of solitary wave packets gathers an increasing amount of the 
correlation in small xenochrysts immersed in uncorrelated low-amplitude 
fluctuations. Up to this point, the process is very similar to the self- 
focusing scenario presented in the spin- and NLS-systems. 
(iii) Rapid divergence: However, despite the fact that for a long time the 
system appears to reach a statistically stationary state, in the end it 
is clear that no equilibrium is attained and the local amplitude rapidly 
diverges. At m 2960, two peaks merge at n ps 420 creating an all- 
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time high of the amplitude that apparently exceeds a certain threshold 
locally. This highest peak now starts to grow rapidly so that the iteration 
is derailed within a few time steps. The features of this rapidly growing 
'monster'-solution will be described in the next section. 

Unlike the modulus square norm of the continuous KdV equation, J2 u i n ) 2 
is not conserved in the whole process. Fig. 9a shows J2( u ( n ) + v(n)) 2 and 
— J2{ u i n ) — v{n)) 2 as a function of the time m. XXw(n) + v{n)) 2 equals 
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Fig. 9. XX n ( n ) + v(n)) 2 and — J2(u(n) — v{n)) 2 of the leapfrog-discretized KdV 
system as a function of time m for the k = 27r/3-wave (a) and for the white noise 
initial condition (b). 

AJ2u(n)v(n) initially while — J2( u ( n ) —v{n)) 2 starts at zero; the sum of both 
quantities ~< u{n)v{n) > is conserved. In the 'reversible' range m < 400 (i), 
both quantities increase and recur to their initial value periodically. As the 
merging of peaks starts (ii), they settle to almost constant values undergoing 
only a very slow increase. This behavior resembles Fig. 4 up to the final blow- 
up (iii) where both quantities diverge rapidly. 

Initial condition with a low amount of < uv > lead to solitary solutions that 
do not reach the threshold for the blow up. Their growth ends when a few of 
them have absorbed this conserved quantity and move with the same speed. 
This state however is also unstable because of an instability to be described 
in the next section. 

3.4-2 Uncorrelated white noise initial conditions 

Uncorrelated low-amplitude white noise initial conditions < uv >= lead a 
creeping nonlinear process that suddenly ends up in the same sort of local rapid 
divergence of the amplitude. The time elapsing until the system blows up is 
inversely proportional to the square of the noise amplitude. For random white 
noise initial conditions with the same amplitude it is Poisson-distributed. 
Fig. 7b shows the spatiotemporal pattern of locations where the amplitude 
slightly exceeds the noise level. Unlike the correlated case, there is no spa- 
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tiotemporal pattern of merging peaks. A change in the amplitude profile of 
the noise background is hardly detectable even shortly before the blow-up 
occurs. Again three phases of the dynamical behavior can be distinguished: 

(i) Regular behavior: For about 32000 time steps, the amplitudes are at 
the level of the noise imposed by the initial conditions. During this time, 
physical structures such as solitons can be simulated reliably when they 
are imposed by the initial conditions. 

(ii) Creeping focusing process A weak non-moving maximum emerges at 
m 33500, n ps 550 and grows slowly. Fig. 10 shows the low-pass filtered 
amplitudes of u(n) (a) and v(n) (b) of this structure at m — 33500, m = 







(a) 




(b) 




ujn) 


0.04 






0.04 


























m=33500y 
\ . 












m=34000y 




















: y;\m=33500 


0.04 










: V: m=M000 














m=3i300\ 




-0.04 


■/ m=34300 



600 n 700 




500 600 n 700 

(d) 



ujn) 




m=4300~-J 
m=4500'j 





Fig. 10. Profile of the peak evolving out of white noise (Fig.7b). Low-pass filtered 
amplitudes u(n) (a) and v(n) (b) at m = 33500, m = 34000, m = 34300. (c): u(n) 
at m = 34374 two time steps before the iteration breaks down. The dotted line is 
the analytical solution, (d): u(n) at m = 4100, m = 4300 and m = 4500 for smooth 
initial conditions vo(n) = 0, uo(n) = 0.01/cosh 2 (n — 510) is very similar to the 
simulation (a). 

34000, m = 34300. The low pass filtered amplitude changes significantly 
with each time step, but very little with two subsequent time steps. Fig. 11 
shows the corresponding spectral density, 
(iii) Rapid divergence After a slow growing process over about 1000 time 
steps this structure suddenly starts to grow rapidly and derails the iter- 
ation. Fig. 10c shows u(n) for this solution at m — 34374. Fig. 11 shows 
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the spectral density of < uv > at the time steps of Fig. 10. This solution 
(a) (b) 



FT(<u(n)v(n+r)>) 



k 7i/3 



FT(<u(n)v(n+r)>) 




k 7t/3 



Fig. 11. Spectral density (defined as the Fourier-transform of the correlation 
< u{n)(v(n + r) +v(n — r)) + v{n)(u{n + r) +u(n — r)) >, r < 100 for the simulation 
of Fig. 7b at the time steps of Fig. 10. (a) shows the slowly growing solution that 
leads to the monster, (b): The correlation for the smooth initial conditions of the 
simulation lOd. 

(which we call the 'monster' solution because its spatial structure resem- 
bles the Loch Ness monster with several undulations of its tail sticking 
out of the water) appears to be the systems canonical trajectory towards 
infinite amplitudes. 

Monsters are strongly localized: Left to the highest negative amplitude, the 
head of the monster, the amplitudes are close to zero. Right to the head, the 
amplitude of the zig-zag tail decreases rapidly. The monster moves to the left 
by one lattice site with each time step. There are also monsters that move to 
the right; their shape is related to left-moving monsters as u(n) — > — u{— n). 
Most significantly, the amplitude of the monster is squared by every time step 
so that the monster grows as exp(exp(m)) . Fig. 7b shows a delta-shaped broad- 
ening zone of high amplitudes near the blow-up since the tail is growing while 
the head moves towards small n. 

This solution can be calculated analytically. For a state u m (n) = Aa{n) with 
a high amplitude A, the linear terms of equation (3) may be neglected. We 
assume that the structure grows quadratically and moves to the right as 
u m+i{n) = 1rA 2 a(n — 1). Setting a(n) = for n > and for negative 
odd n, we get the solutions a(0) = 1, a(— 2) = (1 — \/5)/2, a(—2n — 2) = 
(1 - ^1 + 4a(-2n) 2 )/2 as the solution of a(n - 1) = a(n - l) 2 - a{n + l) 2 . 
The right moving solution is obtained by setting a(n) — > —a(—n). The dotted 
line in Fig. 10c shows this solution where the monster's head is submerged. 
Most significantly, the conserved quantities < u >, < v > and < uv > are zero 
for this solution; the monster needs no external food source for its growth. On 
the other side, it rapidly produces high amounts of < u 2 + v 2 > (Fig. 9b). 
A blow-up with these characteristic features follows from various initial con- 
ditions. In the previous section we have described its evolution out of solitary 
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solutions that emerge from a Benjamin- Feir type of instability and exceed a 
certain threshold after the merging process. It can also grow out of the noise 
between these solitary solutions or out of weak white noise through a much 
more dramatic type of instability. 



4 Statistical analysis of the final state 

In this section we will give a detailed interpretation of these results. The merg- 
ing of peaks corresponds to an Arnold-diffusion process in phase space that 
transfers the trajectory from the initial critical torus to less distinctive parts 
of the shell of constant energy < 7i >= E and magnetization < M. >= M. 
The numerical findings indicate that the characteristics of the final solution 
are determined by the values of the integrals of motion; i.e. that the system 
reaches a thermodynamic equilibrium. One can therefore establish thermo- 
dynamic connections between macroscopic observables of the final solution 
and the integrals of motion. We will compute the equilibrium statistics of the 
system and compare the results to our numerical findings. 

4-1 The partition function 

4-1.1 Low-temperature approximation 

The numerical findings suggest that for low energies the spins point to small 
regions near to the poles (Fig. la) and avoid the region nearer to the equator. 
With a n = ±1 one may approximate spins near to the north- or south pole as 

S nz ^o- n {l- l -(Sl x + Sl y )) (4) 

Low amplitude fluctuations near the north pole are represented by a n = 1 and 
small values of S nx / y . The coherent structures with spins near to the south 
pole correspond to o~ n — — 1. This matching height of all peaks is the main 
technical advantage of the spin chain. Assuming that a n a n+ i = 1 holds for 
almost all n (i.e. the number of domain walls is small), one may approximate 
J2 < 7 n°~n+iSnx/ y ~ En ^nx/y ^ an s P ms are close to the poles, the approximate 
Hamiltonian 

T~teff = En \{Snx + S ny) + J i S nx + S ny) 

—J(S nx S n+ i x -\-S ny S n+ i y ) (5) 
— Ja n a n+ i 
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represents a chain of coupled harmonic oscillators (S nx , S ny ) and a chain of 
Ising spins a n . The approximation holds if the energy is low and the magneti- 
zation is close to the maximum, i.e. J2(l—S nz )/N <gC 1 . It neglects the coupling 
between the oscillators S nx / y and the Ising spin a n so that the Hamiltonian 
splits up into a spin-wave Hamiltonian H w (S nx , S ny ) and an Ising Hamiltonian 
Tti(a n ). H w contains the lowest order terms of S nx / y and neglects anharmonic 
energy contributions of neighboring spins that point to the same hemisphere. 
Tii accounts for the coupling between up and down spins. In terms of the 
stereographic projection, 7i w contains the terms that prevail for |0| <C 1 while 
Hi allows for the contributions for \<p\ 3> 1. The nonlinearity is only reflected 
in the Ising magnet. 

The magnetization as a second integral of motion may be approximated as 
M eff = M T (a) + M W (S X , S y ) = J2°n-Y: \{S 2 nx + Sl y ) (6) 

n n " 

if most of the spins point to the north pole. Again this approximation neglects 
higher order terms in S nx / y and contributions cr n S 2 x ^ y with a n = —1. 

4-1.2 Grandcanonical partition function 

The phase space surface on which M. and H are constant can be computed 
most easily using the grand partition function 

7 ) = J e'^ff-^ffUT (7) 

with two parameters j3 and 7. Unlike the canonical ensemble [19], the grand 
canonical ensemble reflects the second integral of motion by the parameter 7 
that controlls the system's magnetization. (3 is the inverse temperature while 
7 is an equivalent of a magnetic field or chemical potential. The exponent in 
(7) is a sum H. e ff — 1-M-eff — (7~Lw ~ iM-w) + ~ iM-i) of a spin-wave 
contribution 

H w — "yM w = ^2 7;{(Snx — S n+ i x ) 2 + (S n+ly — S n+ i y ) 2 ) H 7T~^~(Snx + Sl y )(&) 

n 1 & 

that depends only on S nx , S ny and an Ising contribution 

TCj — 'jAii = J — 0~n°~n+l) — 7 51 °~n (9) 

n 

that depends only on a n . The partition function (7) of the whole system is the 
product y = y w yi where y w is obtained by an integration over the variables 
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S n xi S ny while yj is a sum of the configurations of the Ising spins a n . The grand 
partition function of the N Ising spins 

Vl = (coshfrP) + fi) N (10) 



with the abbreviation /x = sinh 2 ^^) + e~ 4J/3 is just the canonical partition 
function of an Ising magnet in an external field 7. The linear dynamics of S nx 
and S ny has the symplectic structure S nx / y = ±&H w /dS ny / x so that a phase 
space volume element may be approximated as dT = Y\dS xn dS yn . For /3 ^> 1, 
the grand partition function of the spin waves can be obtained by integrating 
of dT — n dS xn dS yn from minus to plus infinity. Using the abbreviation A = 

\l 2 + + ^ e ^ auss i an integrals yield 

yUP,i) = (^) N (ii) 



4-2 Thermodynamic relations 



4-2.1 Energy and magnetization 

The thermodynamic properties of the equilibrium state may be derived from 
the grand partition function lny((3, 7) = lny w + Inyj. The parameters /3, 7 and 
the conserved quantities < H >= E, < M >= M are connected by 



M M„. • M, 

= -pHHv^) + Hvi)) (12) 

= A^(-^ + =^) 



E = E w + Ei 

= (i*-i)(Mi/J + My/)) (13) 

/3 V A' 005/1(7/3)^+^2 



with A = y 4J(1 + 7) + (1 + 7) 2 . The partition function is valid for small 
energies per lattice site E/N C 1 and for small mean deviations 1 — M/N <ti 1 
of the spins from the north pole. Low amplitude initial conditions without huge 
deviations from the north pole correspond to magnetizations in the interval 
1 - E/N < M/N < 1 - E/(N(1 + 4J)). The lower bound corresponds to a 
spatially homogeneous initial condition while a wave with k = n defines the 
upper bound. 

Mi, Ei, M w , E w are the physically most interesting quantities: 
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Mi is the magnetization of the Ising system and measures the total extent 
of the coherent structures. At its maximum Mj = N, all spins point up 
while smaller values N > Mj > M indicate the existence of coherent 
structures where the spins point down 

Ej is the positive coupling energy of the domain boundaries and determines 

the number of coherent structures 
M w is the negative magnetization of small fluctuations 

E w is the positive energy of small fluctuations and comprises a coupling term 
and an anisotropic term 

These quantities can be found by computing (3 and 7 as functions of E and 
M and then plugging (3 and 7 in the expressions for Ei, Mi, E W) M w . (12), 
(13) can be solved analytically for the low energy case E/N <C 1. The solution 
(Fig. 12) is qualitatively different in the 'oversaturated phase' with M < M eq 
and in the 'overheated phase' M eq < M < M n with M eq = N — E/y/T+JJ: 



Fig. 12. The temperature (3 1 {M W ,E W ) (a) and the chemical potential 7(M [t) , E w ) 
(b) as obtained from the equations (12) and (13). The energy is fixed as E w = 0.1 
in (c). 
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4-2.2 Oversaturated phase M < M eq : 

(i) The temperature ft^ 1 ~ E/N <C 1 is almost independent of M 

(ii) The chemical potential 7 pa Ee~ 2JN ^ E / ^2N(M eq — M) pa is exponen- 
tially small unless the magnetization comes very close to the transition 
point where this approximation breaks down. 

(iii) Almost the total energy E w pa E is absorbed by the fluctuations while 
the surface energy E I ~ e -2JN/E ~ q is exponentially small. 



(iv) The fluctuations share M w pa —E/{yl + 4JN) of the magnetization is 
independent of the total magnetization. The remainder of M — N is ab- 
sorbed by the spins that are flipped down. The share of spins that is 
flipped down is at most of the order of the energy. 



4.2.3 Overheated phase M eq < M < M n : 

(i) The temperature 

! = (E + M-N)(AJ(M-N) + E + (M-N)) 
P (4J(M - N) + 2(E + (M - N)))N 1 ' 

(ii) and the chemical potential 

(E + M-Nf 

1 AJ(M — N) 2 + 2(E + M — N)(M — N) v ; 

both have a singularity at M s = N — E/(2J + 1). The temperature is 
positive between M w = M and the singularity because the number of 
accessible states grows with the energy in this range. Beyond the singu- 
larity more energy leads to a decreasing number of states so that the 
temperature is negative. 

(iii) The spin waves E w pa E again absorb the bulk of the energy while 
Ei ~ e ~ 2JN / E ~ 0. 

(iv) The Ising-magnetization is near to its maximum M/ pa iV independently 
of M. Fluctuations contribute the magnetization M w pa —E/X 



4-2.4 Transition at M t 



eq 



The solution for Mj explains the transition behavior of Fig. 5 (section 3.2.1) 
and the emergence of coherent structures quantitatively. While the oversat- 



25 



urated phase corresponds to long spin-wave initial conditions, the thermo- 
dynamic equilibrium state is characterized by coherent structures, i.e. spins 
pointing to the south pole: Below the threshold M eg , the Ising-magnetization 
Mi deviates significantly from its maximum Mj = N and the number of spins 
that point down increases linearly with M eq — M. Above the transition the 
Ising magnetization deviates very little from its maximum Mj = N, so there 
are no coherent structures. 

In both phases almost all energy is absorbed by low amplitude fluctuations. 
The surface energy Ej is exponentially small, so that the spins form a very 
small number of domains. The higher number of domains obtained numeri- 
cally indicates that the system does not thermalize completely on reasonable 
time scales. 

As the spins interact only pairwise with a short range in one dimension, the 
transition between the two phases is of diffuse type and not a genuine phase 
transition. 7 and Mj are analytic functions. The slope of 7 increases rapidly 
within a small interval ~ e -4JN/E ^ go ^at the transition approaches a 
phase transition as the energy goes to zero. 



4-3 The entropy 

4-3.1 The shape of the entropy function 

The thermodynamic reasons for coherent structures are best described in 
terms of the systems entropy. The conserved quantities M and E are known 
from the initial conditions rather than the arguments f3, 7 of the grand par- 
tition function. Consequently, the entropy as a function of M and E is the 
appropriate thermodynamic potential. The entropy follows from the grand 
partition function by two Legendre transformations 



Both the spin-waves and the Ising system contribute to the entropy. The two 
systems can get different shares Ej, Mj and E w , M w of the two conserved 
quantities E and M. The entropy of the Ising has the form ~ —Ejln(Ei/N). 
The spin- wave entropy per lattice site is given by S w /N = Infl where the total 
number of accessible microstates is Q N with 



S = ln{y) + P(E - 7M) 



(16) 



n = 



(E W + (1 + 4J)M W )(E W + M W ) 
NM W 



(17) 



26 



50 the entropy of small fluctuations depends on the energy as ~ Nln(E w /N). 
Both systems are coupled thermally, so they have matching temperatures fl^ 1 . 
Using (3 = H we reestablish the fact that the Ising energy Ej/N ~ e ~ 2Jf3 is 
exponentially small compared to the energy of the fluctuations E w /N ~ f3~ l 
for low energies. The resulting Ising entropy is again exponentially small 

51 ~ e~ 2JlS compared to the spin-waves contribution S w ~ — lnj3 + const. 
Consequently, the main part of the entropy arises from the degrees of free- 
dom of waves with small amplitudes while the Ising system only provides an 
almost constant contribution. The Ising system can absorb some of the sys- 
tems magnetization without changing its energy and entropy significantly. By 
doing that, the fluctuations share of magnetization can also change allowing 
the fluctuations to maximize their entropy. The maximum of S w = Nlnfl as 
a function of E w and M w is approximately the total entropy maximum. 

Fig. 13a shows the number of states per lattice site f2 of the fluctuations as a 

(a) (b) 




Fig. 13. $1 for the spin chain as a function of E w and M w : (a) focusing equation 
with J = 0.4. The focussing process corresponds to the arrow on the left slope, 
the arrow on the crest sketches the merging of peaks. The arrow on the right slope 
represents the destruction of coherent structures by short- wavelength fluctuations, 
(b) defocusing case for J = 0.1 and J = 0.4. The line M n crosses E w = for 
4 J - 1 = 0. 

function of M w and E w . It has the shape of a crest ascending towards higher 
energies. All possible small-amplitude states corresponding to positive Vt are 
in a triangular region limited by two highly ordered solutions (M = — E w and 
M w = —E W /(AJ + 1)) and by the systems total energy E w = E. 
The area between the lines Mq and M eq again represents the oversaturated 
phase. The rim Mq = —E w with Q = corresponds to a monochromatic wave 
with k = 0. Long wavelength solutions (which are also representative for con- 
tinuous systems) are located at the slope near to this line. Such fluctuations 
have a low ratio Ej (N — M). Solutions that include high peaks may have even 
smaller values M < N + M . However, their inert down-magnetized domains 
have little influence on the thermodynamics and these states are similar to 
the ones at M . 

The overheated phase is located between the lines M eq and M w . The rim 
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M n = —E W /(4J + 1) with Q = represents a wave k = n at the boundary 
of the Brillouin zone. This wave has the highest ratio E/(N — M) of all solu- 
tions. 

The systems total energy E w = E gives a third boundary of the accessible 
states. The absolute maximum of the entropy is located on this boundary at 
M w = -E/y/AJTT. 



4-3.2 Coherent structures 

The formation of peaks (i.e. down-magnetized domains) can be understood 
as the maximization of the entropy under the restriction of the conserved 
quantities. Formation and merging or destruction of coherent structures is 
represented by paths from the slopes to the crest and to the entropy max- 
imum typifying the phenomena that have been observed numerically. While 
these are nonequilibrium processes since M w and E w are changing, equation 
(17) gives the equlibrium entropy for a system thermalizing at particular con- 
stant values of M w , E w . 



(i) Formation of coherent structures in the oversaturated phase: 

For M < N + M eq (or M w < M eq , e.g. for spatially homogeneous initial 
conditions or long waves), the system can increase M w < and decrease 
Mj > by flipping spins from the north to the south. In the entropy 
profile this means that the system is allowed to move from the Mo-side in 
the direction towards M eq (along the arrow at the left slope in Fig. 13a). 
This leads to an increase of the spin-wave entropy S w for initial condition 
on the M -side of the slope. This process stops when the crest is reached 
so that the ideal amount of magnetization M eq is allocated to the spin 
waves. For long-wavelength initial conditions, the formation of coherent 
structures allows the exploitation of short-wavelength degrees of freedom 
to increase the entropy. 

An additional increase of the fluctuations entropy may be reached by 
transferring energy from Ej to E w . This happens when merging down- 
magnetized domains reduce the domain- wall energy contributing to Ej. 
In Fig. 13a we can identify the route along the crest M eq with this process. 
Finally, E w absorbs almost all energy at the summit leaving little energy 
Ej <C E w for domain boundaries. The resulting magnetization of the 
Ising magnet is 

E 

Mj = M+ - (18) 

(ii) Destruction of coherent structures in the overheated phase: For 

M > N + M eq (or > M w > M eq , e.g. a k — 7r-spin wave as initial 
condition), flipping spins down is impossible because this would decrease 



28 



the entropy. The opposite movement starting from the M v slope towards 
the crest M eq is only possible if some spins are already flipped down so 
that they may be flipped up now. This type of thermalization process 
occurs for initial conditions of spin-down xenochrysts immersed in short- 
wavelength low amplitude fluctuations. This process ends if either all 
spins point up [Mj = N) or if the ideal amount of magnetization M eq is 
allocated to the spin waves. The crest of the entropy may be approached 
from the M^-side by melting existing coherent structures away (arrow at 
the right slope in Fig. 3c). 

In Fig. 5 compares numerical and analytical results for the number of down- 
spins as a function of the total energy, while the total magnetization is fixed. 
For low energies (long wave initial conditions), a relatively big number of spins 
points down so that Mj is smaller. For higher energies (smaller wavelengths), 
the number of down spins decreases and reaches zero at the transition point. 
In Fig. 13a, these initial conditions correspond to energies on a line M = const 
connecting points on the lines Mq and M n . The threshold of Fig.5 corresponds 
to the intersection point of the line M = const and the crest M eq (E). This 
threshold is obtained for any path that crosses M eq . Below the threshold, 
the spin-wave entropy can be maximized by flipping spins down according to 
equation (18). Above the threshold only an exponentially small share of spins 
point down. Again, the transition is analytic but very sharp. 
We conclude that the thermalization of energy under the constraint of the sec- 
ond integral of motion produces high-amplitude peaks emerging from an irreg- 
ular low-amplitude background. The formation of coherent structures allows 
the system to increase its entropy of low-amplitude fluctuations by allocating 
the right amount of magnetization to the spin-waves. While the total magneti- 
zation M = M w + Mj is constant, the spin-wave part M w = — ^2{S^ n + Sy n ) /2 
can turn over magnetization to the Ising part Mj = J2 a n that can attain any 
value in the interval M < Mj < N. This enhances the entropy of the spin 
waves while the Ising entropy is negligible. The entropy increases with E w 
so that almost all energy is allocated to the fluctuations. This explains the 
merging process of the peaks where energy is transfered from Ej to E w . 

4-4 Power spectrum and particle conservation 

Spin waves with a wavenumber k contribute the power 

1 

< Uk >_ (3{2J(l-cosk) + l + 1 ) 

to the Hamiltonian E = J2k n k^k with n& = Sk x S~kx + Sk y S-k y an d u k = 
2J(1 — cosk) + I. The 'particle number' X)rifc is related to the magnetization 
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(19) 



M w . Fig.l4a compares the Rayleigh- Jeans distribution < n k >= T/(j + uj k ) 
to numerical simulations. 




Fig. 14. Spatial power spectra of the spin chain averaged over 100000 time steps after 
a previous integration over 100000 time steps and the corresponding Rayleigh-Jeans 
distributions (19). (a) was obtained for homogeneous initial conditions, (b) follows 
from a noisy short-wave (k « 2-7r/3) initial condition related to negative (3 and 7. 



i) In the overheated phase below the transition, the distribution is independent 
of M since 7 is exponentially small. The energy n k LO k is distributed equally 
over the /c-space while small wavenumbers have the highest power n k . The 
power spectrum Fig. 14a is almost unchanged throughout the oversaturated 
phase and the fluctuations energy per particle is constantly y/AJ + 1. The 
systems surplus of particles condenses at the south-pole state with low ener- 
gies per particle, i.e. some spins are flipped to the south pole. The Rayleigh- 
Jeans distribution is attained during the merging process starting from the 
peak-like spectrum of the initial monochromatic wave. The fluctuations de- 
termine the equilibrium ratio of anisotropic energy and coupling energy 
(Fig.4a) as E a /Ej = (VUTT - l)" 1 . Fig.4b compares this formula to the 
results of simulations with various coupling parameters J. Deviations are 
due to the fact that the system does not reach the perfect equilibrium af- 
ter the integration. Some additional domain walls lead to slightly increased 
values of Ej. 

i) Above the threshold, 7 strongly depends on M and E and the power spec- 
trum is deformed. The temperature becomes negative in the strongly 'over- 
heated' domain since the entropy as a function of E decreases for M > M^. 
In this range, most of the energy is due to short waves (see Fig. 14b) with 
/T 1 « -0.2, 7 « -2.7). 
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4-5 Defocusing equation 



Thermodynamics of the 'defocusing' case (Fig.3d) with a negative anisotropy 
T~L — Z)n^(l — S„S n+ i) — (1 — S 2 n )/2 is slightly more complex. We obtain 
= (E + (4J- 1)M W )(E - M w )/M w where E is negative. We distinguish two 
cases: 

(i) For weak coupling 4 J— 1 < 0, f2 has got a maximum at M eq = E/^f—AJ + 1 
(see f2-shell for J = 0.1 in Fig. 13b). Q is zero for Mq = E < and for 
M,r = — Ej (4J — 1) < 0. Differently from the system with a positive sign of 
the nonlinearity in (1), now M n < M , i.e. short wavelength spin-waves have 
the lowest magnetization for a given energy. Starting from M n , the system 
can approach M eq and increase its entropy by flipping down single spins. 
In contrast to the process described earlier, it is now favorable to store a 
maximum of energy in the domain boundaries. 

(ii) For strong coupling (4 J — 1 > 0), Q as a function of M decreases in the 
whole interval of accessible values of M w < M or M w < M n and is zero at 
the homogeneous state M w = M = E and for short waves M w = M w (see 
J = 0.4 in Fig. 13b). The entropy of the spin-waves can not be increased 
by decreasing the magnetization. Thermodynamics allows a weak focusing 
process by storing energy in domain walls starting from M n , but we have 
found this process numerically only in the DNLS system. 



4-6 NLS systems versus spin systems 

Thermodynamics supports the equivalence of spins and NLS-systems with 
respect to the formation of coherent structures. Fig. 15 shows the nonlinear 
NLS-energy — 10„| 4 at lattice sites n as a function of the particle number \(f) n \ 2 
and the corresponding anisotropic spin-energy (1 — S n ^jl'2 as a function of the 
conserved negative magnetization 1 — S nz ; the points indicate typical states 
of the oscillators in equilibrium. The spin system has two states with a low 
energy per lattice site: The north pole state is characterized by high energies 
per particle E/(N — M) while the specific energy is low near the south pole. 
The latter state is more fuzzy for NLS systems since the potential energy 
— \(f> n \ A and the particle number \<p n \ 2 per lattice site are unbounded. Any 
contribution proportional to the second integrals may be added to the energies; 
consequently, the ratio of the potential energy and the particle number per 
lattice site is relevant. In the NLS-system, — |0„| 4 /|0 ra | 2 equals the negative 
particle number at the lattice site. Similarly, the corresponding ratio (1 — 
Snz) /(I — Snz) — 2 — (1 — S nz ) decreases linearly with the 'particle number' 1 — 
S nz . The only relevant difference is again the maximum 'particle number' 1— S z 
per lattice site corresponding to the down-spins of coherent structures while 
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(a) 



(b) 




Fig. 15. Sketch of (a) the anisotropic energy (1 — S 2 z )/2 of the spins as a function of 
1 — S n z an d (b) of the nonlinear energy energy — |</> n | 4 as a function of the particle 
number \(j) n \ 2 . The dots indicate the typical equilibrium position: Most lattice sites 
have low particle densities while small domains attain high particle densities. In the 
spin chain, these domains extend over several lattice sites. 

the peak-amplitude of the NLS-system can obtain various values depending on 
the initial conditions. As a result of that, each domain with a high amplitude 
consist of only one lattice site that absorbs the particles. 
Both systems condense some particles in the low-energetic coherent structures 
in order to increase the energy of the remaining particles. In both systems the 
the energy per particle of the condensates decreases linearly with the particle 
density per lattice site. This energy becomes available to the coupling part of 
the Hamiltonian so that the system can explore degrees of freedom of wave-like 
fluctuations. This disordered state can only absorb a finite amount of energy 
per lattice site since the lattice constant limits the shortest wave lengths. 
In contrast, fluctuations on infinitesimal length scales of continuous systems 
absorb all the energy in the coupling term ~ |V0| 2 as a zero-energy peak 
absorbs all particles so that the system blows up in finite time. 



4-7 Particle-nonconserving systems 

The separation into a high- and a low-amplitude state follows from the en- 
tropy maximization under the restriction of the second conserved quantity. 
This separation of the system into small fluctuations and high peaks does 
not occur if the second quantity is not conserved (section 3.3.4). The system 
can produce or annihilate particles to increase the fluctuations entropy and 
thermalizes on the energy shell without further constraints. The Hamiltonian 
ft = En J^n^n - <Pn<P* n +i ~ 0n0«+i) - 1 n | 2 - 1 0„ | 4 (we neglect the small 
symmetry breaking term) may be considered as a sum E = T + V of a cou- 
pling term and a potential term. The sign of u in —u\(p n \ 2 — |0 n | 4 decides if 
there is a potential well or a maximum at = 0. This explains the two types 
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of thermalization found in section 3.3.4: 

uo < The trajectories change very little under the influence of the weak 
symmetry-breaking field for small times. A pattern of high peaks emerges 
initially, but disappears again since the surplus particles are annihilated 
(Fig. 6c). The system finally settles into a state of small fluctuations trapped 
in the potential well. However, some oscillators may escape from the local 
energy minimum and attain high amplitudes for weaker symmetry breaking 
fields. 

uj > The system thermalizes in a state of high-amplitude fluctuations. In 
the discrete NLS-system, the amplitudes continue to grow without bound 
since particles are created Fig.6d. Energy is transferred from the potential 
to the coupling term so that T and \ V\ both grow. This growth stops stops 
finally so that the amplitudes remain finite. 

It is the shape of the potential that leads to low-amplitude fluctuations by 
particle annihilation (i) or to particle creation allowing high-amplitude fluc- 
tuations (ii) by thermalization. Interestingly, particle nonconservation does 
not led to fluctuations with infinite amplitudes in (ii); the particle production 
stops finally. The reason for this is the mismatch of the orders of the potential 
energy V ~ — 10| 4 and the coupling term that grows only quadratically with 
the amplitude. Beyond certain high amplitudes the coupling energy can not 
absorb any more energy that is released by the potential. A further increase 
of the amplitude would distort the Rayleigh- Jeans distribution towards high 
wave numbers and reduce the systems entropy; particle production has to stop 
therefore. In other words, the energy shell does not contain states with infinite 
amplitudes. 

The Landau-Lifshitz equation with oo < leads to small fluctuations since the 
anisotropic energy has quadratic minima at the poles. The potential energy is 
maximal at the north pole for sufficiently strong values uj > 0, so that large 
but finite fluctuations emerge. 

4-8 Leapfrog- discretized Korteweg-de Vries equation 

The intermediate dynamics of the discretized KdV equation resembles the spin 
systems formation of an equilibrium state. The blow-up however is an intrinsic 
nonequilibrium process. We study the phase space volume that is accessible 
to the system during this process. 

4-8.1 Phase space shell 

The leapfrog-discretized KdV-equation (3) is a nonlinear, area preserving map- 
ping in the N variables u(n),v(n). The area preserving property can be seen 
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from the Jacobian 



J 



■ dG(u(n),...) - 
du(l) 

I 



(20) 



where / is the identity matrix and is the zero matrix. G comprises the linear 
and the nonlinear derivative term of (3). The phase space that is accessible to 
the system is restricted by the intergals of motion J2 u m (n)v m (n), J2 u 2m( n ) = 
Hv2 m +i{n) and E«2m+i(™) = E^2m(n). Introducing the variables P m (n) = 
u m (n)+v m (n) and Q m (n) = u m (n)—v m (n), the conserved quantity J2 u ( n ) v ( n ) 
of (3) may be written as 

2 < uv >= 2 £ u(n)v(n)/N = -L £ P(n) 2 - -L £ Q(n) 2 (21) 



as the difference of two unbounded positive terms while the nonconserved 
modulus square norm is 

< u 2 + v 2 >= 5>(n) 2 + v(n) 2 )/N = -L £ P (n) 2 + -L £ Q(n )2 (22) 



Solutions of (3) which correspond to physical solutions have u(n) w u(n) so 
that ~ 0. Roughly speaking, P(n) is associated to the physical modes 

which are solutions of the original partial differential equation while Q(n) is as- 
sociated to spurious computational modes. The integrals < P >= J2P(n)/N 
and < Q >= J2Q(n)/N are linked to the system's k = modes. It is useful 
to consider the phase space only at even (or equivalently only at odd) time 
steps since the sign of < Q > changes with each time step. 
The modulus-square norm < u 2 + v 2 > is almost constant on an intermediate 
time scale 500 < m < 2800 in Fig.9a and until the blow up (m < 34000) in 
Fig. 9b. While the phase space shell defined by constant < u >, < v > and 
< uv > has an infinite volume, the additional constraint of keeping < u 2 +v 2 > 
constant gives a finite phase space shell v{< u 2 + v 2 >). The total phase space 
shell without this constraint is given by Q ~ / u(< u 2 + v 2 >)d < u 2 + v 2 >. 
The map (3) is again area preserving in P(n) and Q(n), and the microcanon- 
ical partition function can be rewritten as 

v{< u 2 + v 2 >) ~ / 5(E P(n) 2 /2 - N < P 2 /2 >)5(E P(n) - N < P >)U dP(n) 
x / 5(E Q(n) 2 /2 — N < Q 2 /2 >)<5(E Q{n) - N < Q >) n dQ(n) 



The integrals over dP{n) and over dQ{n) each measure the intersection of a 
hypersphere with the radius < P 2 /2 > (and < Q 2 /2 >) and a hyperplane that 
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intersects the P(n)-axes at P(n) =< P >, (and Q(n)-axes at Q(n) =< Q >), 
and one finds 



u(< u 2 + v 2 >) ~ (< P 2 > - < P > 2 ) N -\< Q 2 >-<Q > 2 ) N - 1 ^ 
~ ((< u 2 > + < v 2 > - < u > 2 - < v > 2 ) 2 - 4(< uv > - <u >< v >) 2 ) N ~ 1 



This expression grows rapidly with < u 2 + v 2 >. On the other side, it decreases 
with the correlation < uv >. Since there is no upper limit for the modulus- 
square norm, the surface Q of the conserved quantities phase space shell is 
infinite. 



4-8.2 The focusing process for correlated initial conditions 

During the focusing process Fig. 7a the system gathers huge quantities of the 
correlation < uv > and of < u 2 + v 2 > into dense solitary waves so that 
the remaining space is filled with uncorrelated white noise of u and of v. It 
is tempting to interpret this phenomenon along the lines of the formation of 
down-magnetized domains that allowed the spin chain to maximize the entropy 
of small fluctuations. The partition function (23) however gives no reason for 
this interpretation; < uv > may well be distributed equally in space in (24). 
The main problem is that < u 2 + v 2 > is not conserved. The system can 
increase this quantity to increase its entropy. This is what happens in the first 
400 time steps in Fig. 9a: The initial reversible process increases and decreases 
this quantity as the trajectory is close to a homoclinc orbit linked to the phase- 
instability of the initial k = 27r/3-wave. As the trajectory disappears from this 
orbit, < u 2 +v 2 > reaches a plateau above its initial value that increases slowly 
during subsequent mergings of solitons until shortly before the blow-up. The 
solitons are highly correlated wave-packets of the k = 2n/3 carrier wave, so 
they contain both < u 2 + v 2 > and < uv >. The formation of the solitary 
waves allows the system to increase < u 2 + v 2 > thereby increasing v. 



4-8.3 The focusing process for uncorrelated white noise initial conditions 

While nonlinear terms are extremely weak for the initial noise level ~ 0.01 
(Fig. 7), there is an insidious nonlinear process that increases the amplitude 
slowly to the level where the rapid growth can occur. Fig. 10 shows the low 
pass filtered amplitude of u(n) (a) and v(n) (b) long before the blow up occurs 
at this site. Its shape changes only slowly over even (respectively odd) time 
steps. U2m(n), V2m{n) both grow slowly in time, but they deviate substantially 
from each other. Interestingly, the feedback from short waves changes this 
dynamics very little. Fig.lOd shows the dynamics for smooth long wave initial 
conditions u (n) = 0.01/cosh 2 (n — 510), v (n) = 0. The solution's shape 
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is very similar to the solution emerging out of white uncorrelated noise in 
Fig. 10a and it also initiates the 'monster' solution. The amplitudes may be 
described by continuous functions in space and time as U2 m (n) = r(t,x) and 
tJ 2m (n) = s(t,x). Their dynamics is given by two coupled continuous KdV 
equations 

T* -5 t y y QSSy 

(25) 



These equations describe the long-wave dynamics of the leapfrog-system where 
the strength of the physical and the computational mode is of the same order. 
Unlike the KdV equation, spatially homogeneous solutions of the coupled KdV 
equations can be phase unstable. Constant solutions r = tq and s = sq have 
the eigenvalues A 2 = — r s k 2 — (r + so)^ 4 — k 6 where k has been rescaled. They 
are unstable for r s < and for r = 0, s < (or for r < 0, s = 0). The 
first instability is very similar to the Benjamin- Feir modulational instability 
and leads to traveling solitary waves which eventually can grow sufficiently 
large to access the 'monster' solution. The second is new and particularly 
interesting as it can initiate the 'monster' solution in a region where < uv > 
is zero. Its unstable mode has a shorter wavelength than the Benjamin-Feir 
mode (although still long compared to the lattice constant) and moreover 
grows faster. We argue below that this saddle point in the phase space is 
accessed readily by an evolving solution because of an inverse flux of the 
power spectrum of the u, v correlation towards low wavenumbers. Once the 
system comes close enough to this starting point, a localized solution grows 
irreversibly until it reaches the size necessary for the rapidly growing monster 
solution which is a heteroclinic connection to infinity [20]. 



4.8.4 The blow-up of amplitudes 

As (3) is nonintegrable, it is hardly surprising that the trajectory eventually 
separates from any orbit shadowing a solution of the original partial differential 
equation. The trajectory simply can disappear from the original Kolmogorov- 
Arnold-Moser torus by Arnold diffusion and explore regions of the phase space 
that are most likely connected with high amplitudes. The amazing finding is 
that this process can lead to such a rapid and unpredictable divergence of the 
amplitude. This feature is absent in the spin system where the phase space 
itself is compact. In discrete NLS systems, the coupling energy is restricted 
by the lattice constant and the conserved particle number J2 |0| 2 - The first re- 
striction is absent in the twodimensional continuous NLS equation, a canonical 
example of a system with finite-time wave collapses. The fixed energy is the dif- 
ference of two energies / | V(f>\ 2 dV and / 1 1 4 c/V^ that each are not conserved and 
that can attain any value in a half-open interval. The corresponding conserved 
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quantity in the leapfrog scheme is 2^u{n)v{n) = |l] n -P(^) 2 — \ HnQi. n ) 2 
where J2P( n ) 2 an d J2Q( n ) 2 eac h may grow indefinitely. This suggests that a 
blow-up occurs in systems where 

(i) the phase space non-compact 

(ii) the integral of motion constraining the phase space is the difference of two 
positive unbounded quantities 

We conjecture that the blow-up is the generic way of thermalization in such 
systems. The exploration of the phase space shell defined by a constant dif- 
ference of two positive definite energies leads to finite-time singularities since 
both energies can grow unrestrictedly at the same rate. This may happen if 
a solitary structure grows beyond a certain threshold, or more surprisingly in 
a sudden eruption out of low amplitude fluctuations. The analytically known 
monster solution serves as the canonical highway to infinity during the blow- 
up. 

Analogous to the collapse in the two dimensional nonlinear Schrodinger equa- 
tion this might be an inevitable consequence of a condensation process. In 
that context, the spectral energy whose density is approximately Uk^k (where 
uik = k 2 is the frequency or energy of a wave vector k and rik is the parti- 
cle density or waveaction) has a net flux to high wavenumbers. Because both 
&k n k and Uk are approximately conserved, this means that there must be 
a corresponding flux of particles towards low wavenumbers. This flux leads to 
the growth of a condensate as particles are absorbed at k — 0. For the focus- 
ing NLS equation, this condensate is unstable and leads to the formation of 
collapsing filaments. But the instability is very robust and in fact the collapses 
begin before is all concentrated at k = 0. As soon as there is sufficient n& 
near k — 0, the collapses begin. 

Collapses are inevitable because of the finite flux of particle density to long 
waves. In the present context, a similar scenario occurs. Fig. 11a shows the 
spectral density of < uv > of the structure that leads to the collapsing mon- 
ster. It closely agrees to Fig. lib which shows the spectral density for smooth 
initial conditions of Fig.lOd related to the continuous system (25). The simi- 
larity to the collapse in twodimensional NLS systems suggests that a net flow 
of the spectral density of < uv > towards small wave numbers moves the 
system towards a saddle point for the heteroclinic connection to infinity, but 
the driving force of this process is not yet understood. 



5 Conclusions 

The fusing of the peaks that emerge from an initial phase instability is distinc- 
tive for the self focusing in nonintegrable systems. As opposed to integrable 
dispersive nonlinear systems, it leads to the formation of coherent structures 
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where peaks of high amplitude emerge from a disordered low amplitude back- 
ground. Small wavelength radiation emitted by the fusing peaks leads to an 
irreversible transfer of energy to small scales. This can be understood as ho- 
moclinc chaos at the onset of an Arnold diffusion process where the trajectory 
separates from a critical torus to less distinctive parts of the energy shell in 
phase space. 

This process can be interpreted as the thermalization of energy of an ordered 
initial state. The table 5 summarizes these results. We have discussed the 
equilibrium thermodynamics of a low temperature state that is reached after 
a long time for a generic model, the anisotropic Heisenberg spin chain. The 
system reaches a state where most spins contribute to a disordered state near 
the north pole while the coherent structures correspond to a few xenochrysts 
where the spins point to the south pole (first row in the table). A similar 
state is reached in nonintegrable NLS type of equations but it is absent in 
integrable NLS equations (fourth row in the table) and in systems that do not 
contain a second integral of motion in addition to the energy (third row). The 
latter thermalize in a low amplitude state with no coherent structures while 
integrable systems show continuing quasiperiodic motion. 
The physical conclusion of this result is that the focusing process is driven 
by the generation of entropy in a state of small amplitude waves. The link 
between the entropy and the emergence of peaks is the constraint imposed 
by two integrals of motion. Starting from a highly ordered state, the system 
can not simply increase its entropy by reaching its most likely state of small 
amplitude fluctuations since it is restricted by a second conservation law. To 
reach the entropy maximum, the system has to allocate the right amount of 
the second conserved quantity to low amplitude fluctuations that absorb most 
of the energy. Whenever there is a surplus of this second conserved quantity, 
this maximum can be reached by gathering the surplus of the second con- 
served quantity in the sites of the coherent structures (first row of the table). 
The system cannot thermalize completely and no coherent structures emerge 
if there is a lack of the second conserved quantity (second row) . This leads to 
the phase transition-like dependence of the amount of coherent structures on 
the energy in Fig. 7. 

Speaking in terms of equation (19), an initial distribution of particles n k will 
not be able to reach the Rayleigh- Jeans distribution while obeying both the 
conservation of energy J2 n k^k and the particle number J2 n k- But the restric- 
tion of particle conservation may be circumvented by gathering low-energy 
particles in small domains and transferring their energy to the remaining par- 
ticles to increase the overall entropy. In that sense, the formation of coher- 
ent structures is reminiscent of the condensation of droplets in oversaturated 
steam where the entropy is maximized under the restriction of particle con- 
servation. 

The conserved quantities of the leapfrog-discretized KdV equation are insuffi- 
cient to ensure a thermalization in a well-defined state. This system contains 
new degrees of freedom since the amplitudes at even and at odd times may not 
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Integrals of motion 


Initial conditions 


Phenomena 


Path to maximum 
entropy 


1. Hamiltonian Ti 

2. M (magnetization) 
or J\f (particle number) 
due to rotational 
symmetry 

(section 2.1.3) 


long waves with 
low amplitudes 
low ratio 

< n > / < N > 


Formation of coherent 
structures and low- 
amplitude fluctuations 
(Figs.l, 2, 
3a,b, 6a) 


transfer of < N > into 
coherent structures 
optimum ratio 
of <n> I <M > 
in fluctuations, 
(Figs. 13a, 14a) 


1. Hamiltonian H 

2. M or N due to 
rotational symmetry 
(section 2.1.3) 


short waves with 
small amplitudes and 
mgn < ri > J < N > ; 
coherent structures 
with low 
< n > / < N > 


Destruction of coherent 

structures 

(Fig.3c) 


transfer of < M > from 
coherent structures into 
fluctuations; decrease 
of<H>/<AT> 
in fluctuations 
(Figs. 13a, 14b) 


Hamiltonian 7i 
broken rotational 
symmetry 
(section 3.3.4) 


any 


disordered fluctuations 
(h igs.bc,d) 


production or reduction 

of < N > 

optimum ratio 

of < H > / <AT> 

in fluctuations 


Hamiltonian Ti 
infinite number of 
integrals of motion 
(section 3.3.3) 


any 


quasiperiodic emergence 
of coherent structures 
(Fig.6b) 


no thermalization 


1. energy J2 UV 

2. mass i^ u 7 
(section 2.2.2) 


k = 27r/3- waves with 
nonzero energy 


Formation of coherent 
structures and low- 
amplitude fluctuations 
(Fig. 7a); 

'monster' (Fig. 10c) 


transfer of J2 uv into 
coherent structures, 
heteroclinic connection 
to infinity 
(Fig.8) 


1 . energy uv 

2. mass J2 u, v 
(section 2.2.2) 


white noise with 
zero energy 


sudden local growth 
of the amplitude 
(Figs.7b, 10a,b); 
'monster' (Fig. 10c) 


heteroclinic connection 
to infinity 
(Fig. 10c) 
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be in step. The conserved correlation 2J2u(n)v(n) of these fields is the differ- 
ence of the physical XX M ( n ) +v(n)) 2 /2 and computational J2{ u i n ) —v(n)) 2 /2 
energies that both can grow in an unbounded fashion similar to / \V(j)\ 2 dV and 
/ |0| 4 dV in the collapse of nonlinear Schrodinger systems. An infinite phase 
space volume is accessible on the shells of constant J2u(n)v(n). If ^2u(n)v(n) 
is finite (fifth row of the table), a modulational instability of usual Benjamin- 
Feir type leads to traveling soliton waves which sweep up smaller ones as they 
travel along the circle so that a scenario of merging peaks takes place on in- 
termediate time scales. Finally, the system finds a way to exploit all of the 
phase space when a long-wave instability of the background noise gives rise to 
the rapidly growing 'monster'-solution. The same solution can be initiated by 
a short-wavelength instability for ^u(n)v(n) = (sixth row). This solution 
is the systems canonical way of exploiting the phase space associated with 
infinitely high amplitudes. 
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